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1.  Introduction 


... 


It  has  been  observed  experimentally  (Schubauer  and  Skramstad  1948; 
Ross  et  al.  1970)  that  the  earliest  stage  of  boundary-laye  r transition  on 
smooth  flat  plates  in  low-turbulence  streams  involves  the  rapid  growth  of 
two-dimensional  traveling  waves  known  as  Tollmien -Schlichting  waves. 

With  the  boundary-layer  stability  theory,  only  the  growth  rate  of  these 
plane  waves  can  be  estimated  satisfactorily.  The  necessary  link  between 
the  initial  amplitude  of  these  waves  and  the  ambient  disturbance  level  is  as 
unknown  today  as  it  was  30  years  ago  (Mack  1977).  Tollmien- Schlichting 
waves  can  most  certainly  be  excitable  analytically.  The  nature  of  their 
excitation  has  not  yet  been  defined  but  must  require,  as  does  the  excitation 
of  any  resonant  oscillation,  superposition  of  complementary  solutions  (con- 
ceptually provided  by  stability  theory)  and  of  particular  solutions  for 
boundary- layer  response  to  local  external  disturbances.  Treatments  of 
the  latter  are  incomplete,  and  more  research,  such  as  the  present  theory, 
is  needed  (Reshotko  1976;  Arnal  and  Juillen  1977). 

There  are  three  theories  to  explain  unsteady  infinitesimal  disturbances 
in  a laminar  boundary  layer.  In  one,  it  is  assumed  that  the  surface  is 
vibrated  in  the  streamwise  direction,  producing  uniform  oscillatory  ambient 
flow  (Riley  1975;  Patel  1973).  In  another,  the  effect  of  sound  on  laminar 
boundary  layers  (Illingworth  1958;  Mack  1975)  is  studied,  a coupling  that  is 
fairly  well  understood  but  has  little  effect  in  subsonic  wdnd  tunnels 
(Reshotko  1976).  In  a third  and  in  the  present  study,  the  effects  of  ambient 
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convected  disturbances  are  investigated.  In  these  models,  the  freestream 
turbulence  or  vorticity,  the  most  influential  triggers  of  transition  in  low- 
speed  wind  tunnels,  a e approximated.  Kestin,  Meader,  and  Wang  (1961) 
assumed  that  the  waves  were  oriented  to  travel  purely  in  the  streamwise 
direction.  They  calculated  the  first  departure  from  a quasi-steady  flow. 
Criminale  (1976,  1971)  assumed  that  the  frequency  and  all  wavenumbers 
were  real  and  found  the  effect  of  a point  source  of  normal  velocity  located 
at  the  finite  boundary  layer  edge.  He  concluded  erroneously  that  the 
largest  disturbances  occur  in  the  boundary  layer.  Rogler  and  Reshotko 
(1975)  investigated  a model  in  which  the  imposed  convected  disturbances 
form  a square  vortex  lattice  that  did  not  decay  downstream.  Their  study, 
as  does  this,  indicates  that  maximum  disturbance  velocity  occurs  outside 
the  boundary  layer  and  that  maximum  disturbance  vorticity  occurs  at  the 
surface. 

The  present  theory  advances  the  field  in  several  ways.  Because  it  is 
an  asymptotic  theory  for  very  large  Reynolds  numbers,  closed-form 
expressions  for  maximal  disturbance  velocities  are  achieved,  as  well  as 
explicit  parameter  variations  not  given  by  numerical  calculations  at  a finite 
Reynolds  number.  The  present  formulae  alone  are  simple  enough  that  they 
can  be  combined  later  with  complementary  solutions  from  asymptotic 
stability  theory  to  explain  how  convected  disturbances  generate  Tollmien- 
Schlichtmg  waves.  A fundamental  solution  (for  a harmonic  line  source  at 
a generic  location  upstream)  is  presented,  thus  allowing  the  description  of 


particular  disturbance  fields  such  as  the  square  lattice  of  Rogler  and 
Reshotko  (1975)  by  superposition. 

Recent  studies  of  the  continuous  spectrum  of  the  Orr-Sommerfeld 
equation  (Mack  1976;  Murdock  and  Stewartson  1977;  Jordinson  1971;  Antar 
and  Benek  1978;  Grosch  and  Salwen  1978)  allow  disturbances  of  bounded 
variation  far  from  the  surface.  These  eigenvalue  searches  are  only 
relevant  to  descriptions  of  boundary-layer  response,  but  stability  applica- 
tions were  expected.  In  these  studies,  a finite  boundary-layer  thickness 
was  defined  and  the  Reynolds  number  was  finite. 


[ 


V 
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2.  Freestream  disturbances 

A description  of  the  external  disturbance  field  that  would  exist  in  the 
absence  of  the  flat  plate  helps  to  explain  how  the  theory  for  boundary-layer 
response  to  a singl  wavenumber  component  (Section  3)  will  be  developed 
into  phvsical  variables  (Section  4). 

Only  two-dimensional  external  disturbances  are  considered.  Since 
Tollmien-Schlicht ing  waves  (of  interest  for  later  extension  of  present  theory) 
are  two  dimensional,  it  is  unlikely  that  the  spectrum  of  disturbance  energy 
in  the  third  (transverse)  dimension  is  ot  critical  importance  in  the  initial 
linear  stage  of  Tollmien-Schlichting  wave  generation  Die  spectrum  of 
energy  in  the  normal  and  streamwise  directions  is_  critical  for  this  applica- 
tion, but  there  is  not  a standard  spectrum  in  wind  tunnels.  The  present 
theory  can  be  applied  to  turbulent  streams  by  integrating  the  single- 
frequency theory  result  over  an  appropriate  spectrum. 

In  the  present  theory,  the  external  stream  is  disturbed  only  by  a single 
harmonic  line  disturbance  with  a laminar  wake  (figure  1)  like  a Karman 
vortex  street.  Of  course,  in  the  usual  street  scene,  the  vortices  are  dis- 
crete and  have  interacted  to  form  a double  row.  Linearity  requires  colinear 
vortices  and  permits  continuous  shedding. 

There  are  several  reasons  why  this  disturbance  was  chosen.  It  is  the 
most  fundamental  harmonic  disturbance  from  a line  source,  and  any  two- 
dimensional  disturbance  field  can  conceptually  he  generated  by  use  of  the 
present  result  as  a Green's  function.  This  disturbance  model  closely 
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Figure  1.  Externally  Disturbed  Boundary  Layer 


simulates  the  effects  of  the  vibrating  ribbons  used  in  the  critical 
experiments  (Schubauer  and  Skramstad  1948;  Ross  et  al.  19  '0;  Kachanov, 
Kozlov,  and  Levchenko  1976).  It  enables  a physical  interpretation  of  the 
results  (disturbance  maxima  are  presented  in  Section  4 in  physical  variables, 
rather  than  in  wavenumber  only).  It  also  provides  justification  for  studying 
only  wavenumbers  whose  order  of  magnitude  is  comparable  with  reciprocal 
boundary -layer  thickness . 

If  x denotes  the  distance  downstream  from  the  leading  edge,  and  y the 
normal  distance  from  the  surface,  a disturbance  stream  function  t,  provides 
disturbance  velocity  components  and  A small  cylinder  or  ribbon 

(one  with  a width  that  is  negligible  compared  to  its  distance  from  the  leading 
edge)  is  located  at  Xq  and  v^.  If  at  time  tQ,  the  disturbance  source  sheds  a 
single  vortex  of  circulation  T,  the  vortex  must  consist  of  both  a potential 
vortex  and  a rotational  core,  or  the  vortex  center  will  not  exhibit  solid-body 
rotation  (Rott  19t>4).  The  rotational  core  must  be  convected  with  the  flow  in 
accordance  with  the  vorticity  transport  equation,  and  the  potential  vortex 
must  accompany  it,  on  physical  rather  than  mathematical  grounds.  The 
disturbance  stream  function  in  the  absence  of  the  plate  then  is 


and  Ej  denotes  the  exponential  integral.  If  such  vortices  are  shed 
continuously  and  harmonically,  beginning  at  time  tQ  = 0,  with  frequency  uu 
and  with  peak  circulation  addition  rate  T , then  ijj  is  the  real  part  of 
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The  real  part  of  ip  is  assumed  throughout.  If  T is  introduced  as  the  time 
since  turnon , 
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Only  the  steady-state  velocity  components  are  of  interest.  They  are 
obtained  by  differentiating  0 with  respect  to  y or  x and  taking  the  limit  as 
t ‘ 00 « which  exists.  When  the  Reynolds  number  that  is  based  on  the  dis- 
tance downstream  from  the  source  x - x^  is  very  large,  the  velocity  com- 
ponents are  explicitly 


exp 
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whe  re 
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(b) 


The  error  is  of  order  one  over  the  root  Reynolds  number.  It  can  be  verified 
that  dt'/dy  is  odd  in  y - y and  dv/dx  is  even  in  y - y ^ . Both  decay  tran- 
scendentally  as  the  magnitude  of  y - y^  grows  without  bound. 

It  is  desired  to  examine  the  y-Fourier  transform  of  these  components, 
because  the  theory  for  disturbances  in  the  presence  of  a plate  is  expressed 
as  the  response  to  a single  transverse  wavenumber.  The  velocity  trans- 
forms are  defined  by 


(7) 


and  found  to  be 
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The  argument  of  the  exponential  function  should  have  a real  part  of 
order  one  if  there  is  to  be  significant  disturbance  amplitude.  In  addition, 
it  is  not  limiting  to  assume  that  x - x^  is  comparable  with  the  distance  x 
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from  the  leading  edge  of  a plate,  if  a plate  were  present.  Therefore, 
normalized  frequency  ii  and  transverse  wavenumber  L are  defined  as 
follows. 
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The  parameters  and  L are  assumed  to  be  of  order  one.  L is  real;  is 
positive.  R is  the  large  parameter  in  this  asymptotic  theory.  In  most 
stability  theories,  R is  the  Reynolds  number  that  is  based  on  some 
bounda  ry  - laye  r thickness,  but  here  it  is  the  root  x-Reynolds  number  for 
simplicity.  With  normalized  parameters,  (8)  becomes 


i*IV’UL!  ipHr  -P 

U<*>l  |nl  UR,“  * L 1 


(x  - V 2 2 

(sr  + L ) - iLTl0 


where  "H  q is  y^R/x.  Thus,  the  laminar  wake  behind  an  oscillating  vorticity 
source  grows  parabolically  downstream.  Its  width  is  of  order  x - Xq  over 
the  root  Reynolds  number  based  on  x - Xq.  This  characteristic  thickness 
establishes  the  characteristic  wavenumber  region  of  interest  to  be  that 
where  L is  of  order  one.  Since  the  exponential  function  argument,  in  square 
brackets  in  (10),  has  a negative  real  part  of  order  L^,  no  larger  j L than 
those  of  order  one  need  be  considered  for  this  type  of  disturbance  source. 

Velocity  transforms  for  other  line  disturbance  sources  can  be  found. 
For  example,  if  the  line  source  is  a symmetric  cylinder  that  oscillates  in 


the  streamwise  direction,  it  may  develop  a periodic  drag  force  by  shedding 
a vortex  doublet  sheet.  Velocity  transforms  for  such  a source  also  indi- 
cate exponential  behavior.  However,  if  the  line  source  does  not  shed  any 
vorticity,  but  excites  only  irrotational  disturbances,  then  these  disturbances 
die  out  so  fast  downstream  that  they  are  of  no  interest.  For  example,  the 
stream  'unction  for  a stationary  source-sink  oscillator  can  be  found,  dif- 
ferentiated, and  transformed.  The  y-transforms  of  the  velocity  components 
are  proportional  to  e.xp  [-  (x  - x ) x 1 L R|.  This  bracketed  exponent  is  a 
factor  of  R more  negative  than  that  in  (10),  where  the  disturbance  is  a combi- 
natn  n of  rotational  and  convected  potential  flow. 

In  the  absence  of  a plate,  the  disturbance  velocity  components  are 
recovered  by  taking  the  inverse  transform 
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In  the  presence  of  an  aligned  semi-infinite  flat  plate,  the  disturbance 
velocities  are  more  complicated.  Additional  terms  are  needed,  whose  form 
in  the  freestream  can  be  found  from  the  Oseen  form  of  the  vorticity  trans- 
port equation.  Thus,  outside  the  plate  boundary  layer  the  velocities  are 
given  by 
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where  a fourth  term,  oc  exp(ujy  / U) , has  been  disallowed  on  physical  grounds. 

»>-  u. 

The  terms  containing  C and  C are  for  reflected  and  companion  waves, 

r w 

respectively,  which  are  generated  by  the  presence  of  the  no-slip  plate.  For 
the  description  of  the  inside  of  the  plate  boundary  layer,  the  terms  in  paren- 
theses are  replaced  by  functions  of  boundary-layer  variables  (Section  3). 

Use  of  (12)  for  the  total  external  disturbances  has  several  implications. 
Since  integrals  over  in  the  range  0 <.  x^  < x are  avoided,  it  is  assumed 
that  disturbances  near  x will  not  include  residual  disturbances  created  as 
eigenmodes  in  the  boundary  layer  near  x^  upstream.  Equation  (12)  provides 
the  forced  response,  but  not  the  eigenmodes.  A description  of  eigenmode 
generation  along  the  plate  requires  a subsequent  study. 

Insertion  of  the  plate  on  the  point  disturbance  source  wake  displaces 
the  ambient  flow  and  modifies  the  station  at  which  the  wake  centerline 
enters  the  plate  boundary  layer.  This  correction  for  nonparallel  flow  is 
needed  in  the  freestream,  where  x - x^  is  of  order  x,  but  is  not  needed 
inside  the  boundary  layer,  where  the  flow  is  parallel.  Outside  the  boundary 
layer,  the  ambient  velocity  components  arc  assumed  to  be  (U,  0),  when,  in 
fact,  they  are(U,  U0/2R),  where  0 is  a constant  (^  1 .7208)  associated  with  the 
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Blasius  function.  Therefore,  a source  actually  located  at  x^  and  is  moved 
to  the  apparent  location  x^  and  y^  + Sx/R,  from  whence  the  above  equations 
can  be  applied.  Only  when  the  source  wake  is  too  distant  to  have  met  the 


boundary  layer  at  x does  this  simple  shift  fail. 
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3.  Boundary -layer  disturbances  for  any  external  source 

The  disturbance  stream  function  caused  by  the  contact  of  a wave  of 
given  wavenumber  with  a boundary-layer  edge  is  assumed  to  have  the  form 
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in  the  boundary  layer.  Only  one  value  of  streamwise  wavenumber  k is  use- 
ful for  given  1 and  io.  That  value  is  k s:  oj/U,  if  x - x^  > 0.  Substitution 
into  the  Oseen  form  of  the  vorticity  transport  equation  indicates  that  this 
is  one  of  two  possible  values.  The  other  pertains  when  x - Xq  *■  0 (down- 

stream source  with  upstream  influence,  a problem  not  addressed  here). 

A normalized  streamwise  wavenumber  component  K is  defined: 
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It  is  well  known  that  a leading  approximation  for  o(  ) that  is  uniformly 
valid  for  values  of  P,  up  to  order  one  can  be  obtained  for  this  parameter 
domain  by  solving  the  Or r-Sommerfeld  equation. 


(h  - r)  ("nn  - K'")  - hnn 0 taw 


-K  o, p +■  K 
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(15) 


where  f is  the  Blasius  function,  defined  as  mean-stream  function  over 
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In  fact,  solutions  of  the  Orr-Sommerfeld  equation  are  valid  descriptions 
for  any  T)  less  than  order  InGR.  Fortunately,  the  largest  T)  that  must  be  con- 
sidered that  is  not  in  the  freestream  is  only  of  order  (lnGR)^2.  Boundary 
conditions  are  o(0)  o-^(O)  = 0,  and  as  ri  — ~ 00 , <p  exp  (iL1!)  plus  additional 

waves  of  bounded  variation. 

The  main  advantage  of  the  present  asymptotic  theory  for  large  R is 
that  not  all  terms  in  (15)  are  needed  everywhere.  Thus,  the  following  flow 
layers  are  described  separately:  a wall  layer,  where  T)  is  of  order  (GR)  , 
the  main  boundary  layer  where  T|  is  of  order  one;  an  edge  layer,  where  11 
is  of  order  (InGR)^2;  and  the  freestream.  Only  the  leading  approximation  in 
each  subregion  is  of  interest.  (The  perturbation  term  of  order  R 1 in  (14)  is 
needed  only  in  the  edge  layer,  where  it  provides  a first-order  estimate  of  the 
contents  of  the  first  parentheses  in  (15).  Elsewhere,  K as  G.  ) 

3.1  Wall  layer:  T1=0(GR)  ^ 

As  in  the  asymptotic  stability  theory,  there  is  a thin  layer  adjoining  the 
surface  in  which  the  highest  derivative  is  needed  in  order  to  satisfy  homoge- 
neous conditions  at  the  surface.  To  leading  order  in  R , (15)  reduces  to 


-0-n-n  = <i0R> 


when  1 is  of  order  (OR)"1^2.  The  solution  that  satisfies  o(0)  = 0T(0)  - 0 
and  that  does  not  grow  exponentially  as  1]  approaches  order  one  is 


O C j-1  + (1  - i)  H + exp  [( - 1 + i)  T\  ]j 


The  scale  factor  C will  be  determined  by  matching  with  the  boundary-layer 
function.  From  this  sublayer  profile,  the  transformed  surface  shear  is 
given  by 


iCS2R~ 

o (0)  = - 5 — 

yy  l 

’ ’ X 


The  maximum  disturbance  velocity  component  in  this  layer  occurs  in  the 
streamwise  direction  when  ’(S2R)  ^ = 3.23.  That  value  is 


o 1.069  C (12R)1/2  — 

y.max  x 


The  wall  layer  profile  of  normalized  0 is  shown  in  figure  2.  This  minor 
maximum  is  shown  in  figures  5 and  6 of  (Rogler  and  Reshotko  19  '5)  also. 

3.  2 The  main  boundary  layer:  3]  = 0(11 

When  both  31  and  f__  are  of  order  one,  the  right  member  of  (15)  can  be 

>1  'I 

omitted,  and  the  governing  equation  is  known  as  the  Rayleigh  equation  (with 
phase  speed  equal  to  ambient  speed): 


(f-r  ” 1)  ( " il  <?)  * ^ "R  TJ  ® 


In  order  to  match  the  sublayer  profile,  it  is  required  that,  when  \ is  of 
order  one  and  T)  — 0 , 


C (1  - i) 


i 


V IX1R) 


The  scale  factor  C remains  to  be  found  by  matching  the  calculated  boundary- 
layer  profile  with  that  in  an  edge  layer  to  be  described. 


Numerical  integration  of  the  Rayleigh  equation  indicates  that  the 
velocity  components  both  increase  monotonically  in  the  main  boundary  layer 
with  the  streamwise  component  always  larger  than  the  normal  component. 
Limiting  profiles  of  qT,  (proportional  to  streamwise  velocity  component)  are 

1(1  - fp)_1  - f f (l-f.)'2dT|  U « 1 (22a) 

0 

, 

PT  (0)  ~ I 

Vcosh  ill  Q » 1 (22b) 


These  and  intermediate  S2  profiles  are  shown  in  figure  3.  A second 
approximation  in  the  low-frequency  limit,  for  which  (22a)  is  the  first 
approximation  will  be  used  in  figure  4 and  is  developed  in  Appendix  A. 

At  the  outer  edge  of  the  boundary  layer,  the  Blasius  mean  velocity  and 
shear  profiles  have  an  exponential  asymptotic  form  given  by 
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Figure  3.  Strcamwise  Disturbance  Profile  in  Boundary  Layer 
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where  EI=(T1  - 3)/2  and  TST  refers  to  transcendentally  small  terms,  is 
a Blasius  constant  a:0.  23373.  When  these  coefficient  functions  are  used  in 
(20),  o approaches  a parabolic  cylinder  function  as  T|  — Specifically, 


O ~ D 


(2hV 


( 2 1 1 


2,3 


(25) 


where 


1 Pj  2 + (r.“  - l ) 

8 p,  (C4  - 5h2  + 8)  Pl  - 10 
12  p ^ = (Q4  - 9C2  4-  22)  p2  - 10  pt  + 74 


and  where  D is  introduced  as  a scale  parameter  for  the  outer  edge  of  the 
boundary  layer.  It  must  be  related  to  C,  its  inner  edge  counterpart,  so  that 
disturbances  in  the  wall  layer  can  be  scaled.  D will  be  found  by  matching 
(25)  with  its  corresponding  expansion  obtained  from  the  edge  layer.  The 
ratio  C/D,  obtained  by  numerical  integration  of  (20),  is  plotted  in  figure  4. 
Because  (20)  has  real  coefficients,  C/D  is  real.  Also  shown  are  the  low- 
and  high-frequency  forms  obtained  in  Appendix  A. 

3 . 3 Edge  layer:  T)  -0[ln(GR)]*2 

Equation  (20)  is  inadequate  when  , becomes  logarithmically  large; 
specifically,  when  H ~ H , where 


(20) 


, H 

! r e ' \ j • R » 1 


There,  the  terms  La  the  Rayleigh  equation  have  become  as  small  as  the 
largest  viscous  term  in  the  Or r -Somme r feld  equation.  Ile  and  Ti  + 3 

are  plotted  in  figure  5. 

> 2 

The  sub  region  where  II"  is  in  the  neighborhood  of  11“  is  called  an  edge 

layer.  It  can  be  considered  to  be  simply  a critical  layer  from  asymptotic 

stability  theory  in  the  limit  where  the  phase  speed  has  increased  to  equal 

the  ambient  speed.  The  dimensionless  edge-layer  thickness  here  is  much 

larger  than  the  thickness,  &-  ( R)  , of  an  imbedded  critical  layer  because 

of  the  exponential  tail  of  the  Rlasius  flow  profile.  With  the  present  theory, 

the  edge  - la  ye  r problem  \s  solved  tor  the  tirst  time.  A previous  attempt 

was  made  by  C'.raebel  (1966),  who,  bv  truncating  the  mean  flow  profile  at 

y , concluded  erroneously  that  the  thickness  ot  the  edge  critical  layer 
-1/4 

was  of  order  ( R) 

When  11“  is  of  order  log(.R),  the  Or  r-Somme  rfeld  equation  can  be 
reduced  to 
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The  least  trivial  problem  results  when  the  edge-layer  variable  is  chosen 
to  be  h,  where 


2 2 
H = H + h 
e 


and  h is  of  order  one.  Then,  in  the  edge  layer,  0 is  represented  as  an 


expansion  of  the  form 


0~On(h)  + H"20.(h)  + H'40,(h)  + 
U el  e Z 


where  Oq  and  pj  must  satisfy 


pO,hhhh  + le  (o0,hh  ‘ V " ° 


-h  2 2 

0,  ui  ..  + ie  (O.  I,  - 0. ) - 2h0.  , . , , - 3o_  . . . - ( L -0)0... 

l.hhhh  l,hh  1 ^O.hhhh  O.hhh  0,hh 
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liquation  (29)  is  in  inverse  powers  of  , but  the  true  error  is  even 


larger,  of  order  H * since  the  Orr-Sommerfeld  equation  lias  already  omitted 


a nonparallel  flow  term  containing  the  third  derivative.  The  missing  term 


is  the  VO  term  from  the  vorticity  transport  equation.  This  omission 

yyy  ' 


means  that  the  parallel-flow  assumption  in  this  application  carries  an  error 


of  order  flog(QR)]  ^2.  A numerical  estimate  of  the  size  of  this  error  can 


I 


be  obtained  by  substitution.  Thus,  when  0 is  one  and  the  x-Keynolds  number 

is  the  minimum  critical  value  from  stability  theory  (60,000),  [log(QR)]  * ^ 

is  0.  43.  Disturbance  profile  shapes  in  other  layers  do  not  carry  such  a 

large  error,  but  disturbance  amplitude  levels  do. 

The  lUasius  function  f„.  for  mean  velocity  curvature  is  not  a valid 

estimate  if  ft  is  too  large.  The  mean  curvature  is  needed  as  a coefficient  in 

the  Orr-Sommerfeld  equation.  Goldstein  (1956,  I960)  and  Imai  (1957)  have 

-2  2 

shown  the  error  in  Blasius  functions  to  be  of  order  R log  R when  t is  of 

order  one.  It  seems  reasonable  that  the  relative  error  in  f^p^,  where  fp,, 

2 2 

decays  as  exp  (-11  ),  may  reach  order  one  when  H is  approximately 
->  2 

log  (R  ).  In  the  edge  layer,  however,  11“  is  only  order  log(»  R);  therefore, 
the  relative  error  in  t'^p,  there  should  be  no  more  than  order  roughly  1/R. 

Four  solutions  for  the  fourth-order  equation  (30)  with  variable  coeffi- 

. A D (1)  (2)  (3)  . (2) 

events  are  developed  in  Appendix  13.  these  are  Oq  , Oq  , Oq  , and  Oq. 

The  first  three  are  found  by  recognizing  expansions  about  h The 

fourth  is  derived  from  an  expansion  about  h - -<®.  As  h -*  -00  (boundary- 

layer  side),  the  four  edge-layer  solutions  become 
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where  y is  Euler's  constant  (xO.  577).  In  this  limit,  p^  and  o^3*  grow  too 
rapidly  to  be  separately  useful  in  the  boundary  layer.  These  two  solutions 
asymptotically  satisfy  the  differential  equation 


°0, hhhh 
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instead  of  the  second-order  equation 


(36) 
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which  is  compatible  with  the  bounda  ry- laye r expansions.  Eor  this  reason, 
and  for  compatibility  with  the  boundary-layer  scaling  parameter  D,  the 
four  solutions  are  combined  (weighted)  as  follows: 


~ E |oq"  + (ni  - 2 - 4y ) Oq2)  + (6.  46404  - i 6.  76836)  o^1  ’ j 


+ D y CRH1'  ~3  U)0n 
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Coefficients  D and  E remain  to  be  found  by  matching  with  freestream 
expansions. 

As  h — + oo  (freestream  side),  the  four  edge-layer  solutions  become 
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(39) 


<^2)  - h + TST 


(40) 


(3)  .2  , 

>1  ~~  h 4- 


Oq  ~ h + TST 


(41) 


(2)0  ~ (i/3)  h3  + e2h2  + e + eQ  + TST 


(42) 


where  e^,  e ^ , and  e,  are  constants  that  are  given  in  Appendix  B.  Then,  the 
leading-order  solution  can  be  written 
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Solutions  of  higher  order  were  not  sought:  nevertheless,  from  (31)  and  a 
corresponding  equation  for  the  asymptotic  behavior  as  h-“>  can  be 
established  to  be 


- 6f3  (L2  - a2)  jr  - [i»f3  + 2f2  (L2  - 0 


2>]5? 


+ cubic 


(48) 
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(49) 


3.4  Freestream:  T|  - 7)e  = 0(1) 

Equation  (27)  contains  the  few  terms  needed  in  the  ambient  flow  where 
the  differential  equation  has  constant  coefficients: 


°T]T|T)T)  + (L  ‘ 0 ’ °T)T|  ' " L 0 0 


(50) 


The  solution  desired  is 
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0 ~ e 


iL(T)  - 11  , 
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e+Ce  e+Ce  e 

r w 


(51) 


The  parameters  and  are  more  convenient  than  and  C w used  in 
(12).  The  coefficient  of  the  first  modal  function  [exp  (iLf))]  has  been  nor- 
malized as  indicated  after  (15).  A fourth  modal  function  [ ~ exp(QT))]  has 
been  omitted  in  order  to  satisfy  the  other  of  two  boundary  conditions  at 
namely,  that  q|  not  grow  with  T)  in  the  freestream.  The  second  and  third 
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modes  retained  can  be  interpreted  as  a reflected  wave  of  relative  amplitude 

C'  to  be  determined  and  a potential-flow  response,  which  is  a streamwise 

wave  entirely.  Lack  of  decay  as  Tl  — ® in  the  amplitude  of  the  reflected  wave 

is  acceptable  because  the  exponent  depends  on  wavenumber,  which  is  the 

integration  variable  for  the  inverse  Fourier  transform. 

In  order  that  the  edge-layer  expansion  given  by  (29),  (43),  (48),  and  (49), 

match  (51),  the  edge-layer  expansion  is  written  in  the  freestream  variable 

T|  - T|  by  means  of 
e 


h H 


V Mi) <’  • V 

Upon  substitution,  the  edge-layer  expansion  becomes 
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Only  two  orders  iind  have  been  retained  because  only  the  leading- 

order  solution  in  Appendix  B was  developed,  and  higher-order  solutions 

appear  at  the  next  order.  Equation  (53)  must  match  the  expansion  of  (51) 

for  small  f|  - T,  , which  is 
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L r w. 
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Since  (53)  does  not  contain  terms  containing  (T|  - T]  ) or  (T)  - T|  ) , such 
terms  in  (54)  must  vanish  by  a proper  choice  of  C and  C , thus 


U + iL 
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In  addition, 
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and 
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were  chosen  to  match  powers  of  (T)  - T)  ) and  (T)  - T|  ) . Higher  powers 
(four  through  seven)  also  match  to  leading  order.  The  asymptotic 
variation  of  the  higher-order  functions  dj  and  was  unnecessarily 
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displayed  in  (53)  as  a confirmation  of  the  adequacy  of  the  matching 
procedure.  Substitution  of  (56)  and  (57)  into  (44)  and  (45)  indicates  that 
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and  that  none  of  the  constants  e 7 , e^,  or  e^  is  needed. 

When  coefficients  E and  D are  replaced  in  (38),  in  the  edge  layer  proper, 
where  Oq  ' are  all  of  comparable  magnitude,  the  term  with  E is  theoretically 
larger  than  that  with  D,  by  a factor  of  the  "large"  parameter  H . Thus,  E 
scales  the  edge-layer  disturbances,  and  (C/D)D  scales  disturbances  in  the 
boundary  layer  and  sublayer.  The  leading-order  edge-layer  solution,  from 
(38)  and  Appendix  B is,  therefore,  the  closed-form  expression 

-p-  (6.  46*404  - i 6.  76836)  (1  + ie'h)  + f(2e'h/2)  (60) 

where  f(x)  is  given  by  ( B - 37).  The  magnitude  and  phase  of  o^/E,  which  is 
proportional  to  the  transverse  velocity  component,  and  ^/E,  w^c*’1 
proportional  to  the  streamwise  disturbance  velocity  component,  are  plotted 
in  figures  6a  and  6b,  respectively.  The  maximum  disturbance  velocities 
occur  in  the  freestream,  not  in  the  edge  layer.  Further,  the  vertical  com- 
ponent of  disturbance  velocity,  which  was  numerically  negligible  compared 
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to  the  horizontal  component  in  the  wall  layer  and  main  boundary  layer, 
finally  becomes  comparable  with  the  horizontal  component  and  may  exceed 
the  latter  in  size  in  the  outer  portion  of  the  edge  layer,  depending  on  their 


relative  scales. 

Figure  7 is  a composite  plot  of  the  streamwise  perturbation  velocity 
profile  in  the  three  layers  into  which  the  steady  boundary  layer  is  divided  by 
this  asymptotic  theory,  o-  (which  depends  on  y but  not  t ) is  multiplied  by 
a transform  (which  depends  upon  f but  not  y)  and  then  integrated  over  £ 

to  obtain  a dimensional  velocity.  1 he  integral  of  T di  ser\es  merely  as 
a scale  factor,  one  that  depends  on  the  nature  of  the  external  disturbance. 
The  specific  relationship  between  the  three  sublayer  profiles  depends  only 
on  f.,  R,  Xq,  and  y^.  Experimentally  obtained  profiles  with  a vibrating 
ribbon  (Kachanov,  Kozlov,  and  Kevchenko  19  *o)  and  numerical  results  with 
a convected  lattice  (Rogler  and  Reshotko  1975)  (figures  5 and  6)  both  exhibit 
the  features  of  this  composite  profile,  i.  e.  , a slight  overshoot  in  a thin  sub- 
layer and  very  rapid  growth  with  increasing  distance  trom  the  surface. 


4.  Results 


Disturbances  in  the  presence  of  a semi-infinite  flat  plate,  generated 
as  described  in  Section  2 by  an  external  oscillating  source  of  vorticity  (and 
limited  to  the  forced  response  as  opposed  to  eigenmodes),  are  plotted  for 
two  locations:  the  treestream  and  the  wall-layer  edge. 

In  these  applications,  it  is  necessary  to  adjust  the  disturbance  scale 
differently,  depending  on  which  disturbance  velocity  component  is  desired. 
Specifically,  the  transverse  disturbance  velocity  in  the  presence  of  the 
plate  is 


exp 


x " xo\l  dt 

U / dx 
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since  the  limit  o(^l)  — expfiLT])  plus  bounded  terms  is  required  in  the  free- 
strearn  and  since  dt/dx  is  proportional  to  q(f|)  rather  than  its  derivative. 

The  streamwise  component,  dt/dy.  proportional  to  o-y  is  obtained  from 
(51)  of  Section  3.  Ln  the  freestream,  the  form 
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is  compatible  with  (11)  if  the  streamwise  component  is 
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Equation s (61)  and  (63)  correspond  to  (11),  which  applied  in  the  absence  of 
the  plate. 

4.  1 Freestream  profiles 

Profiles  of  disturbance  velocity  components  in  the  freestream  are 
desired  both  for  comparison  with  recent  experimental  profiles  obtained 
behind  a vibrating  ribbon  (Kachanov,  Kozlov,  and  Levchenko  1976)  and 
because  disturbances  there  are  the  largest.  The  combination  of  (51),  (55), 
and  (61)  through  (63)  resulted  in  the  closed-form  profiles: 
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whe  re 
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and  where  C’j,  S^,  are  integrals  defined  by 
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and  evaluated  by  parameter  differentiation. 


The  closed-form  results  are 


Streamwise  perturbation  profiles  are  compared  in  figure  8 with 
experimental  measurements  obtained  (Kachanov,  Kozlov,  and  Levchenko 
ll*7t,|  behind  a vibrating  ribbon  in  the  presence  of  a plate.  As  described  in 
a Future  at  Virginia  Polytechnic  Institute,  the  experiment  was  performed  at 
a frequence  of  11-1  11/  and  an  ambient  speed  of  532  cm/sec.  The  ribbon 
location  was  not  specified.  (Its  location  is  not  critical,  and  it  has  been 
assumed  to  lie  above  the  leading  edge.)  The  mean  bounda  ry  - laye  r "edge'1 
and  surface  were  located  at  6,  and  S,  respectively.  The  theoretical  y^  - y^ 
(the  minimum  coordinate  plotted)  was  optimized.  The  arbitrary  amplitude 
scaling  is  the  same  at  all  stations. 

Figure  8 indicates  that  the  characteristic  diffusion  width  is  only  half  of 
the  theoretical  value  and  that  the  decay  is  much  more  rapid.  Since  these 
discrepancies  would  also  occur  in  the  absence  of  a plate,  either  (5)  is  wrong, 
which  seems  most  unlikely,  or  the  experimental  frequency  must  have  been 


-uot  5[.a  v EQ.64  ° v-  EXPERIMENT  J 

x = 150  mm  250  mm  350  mm 

Figure  8.  Streamwise  Disturbance  Profiles  Outside  Boundary 
Layer  (w  = 716  rad/sec,  v = 532  cm /sec) 
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twice  the  reported  114  Hz.  If  it  is  the  latter,  quantitative  as  well  as 
qualitative  agreement  could  be  reported. 

4.2  Wall-layer-edge  disturbance  amplitude 
Infinitesimal  external  disturbances  can  only  cause  boundary-layer 
transition  to  the  extent  that  they  excite  boundary-layer  instabilities.  Since 
Tollmien-Schlichting  waves,  i.e.  , two-dimensional  eigenmodes,  travel  at  a 
phase  speed  of  about  one-third  the  ambient  speed,  corresponding  to  a 
critical  layer  (of  asymptotic  stability  theory!  located  near  T|  = 1 , as  much  or 
more  interest  must  be  focused  on  the  forced- response  disturbances  of  rela- 
tively small  amplitudes  near  7|  1 , as  in  the  very  much  larger  amplitudes 

in  the  freestream.  In  this  central  portion  of  the  main  boundary  layer,  the 
st  reamwise  component  is  always  much  larger  in  magnitude  than  the 
transverse  component  and  is  found  (figure  31  for  any  )|  once  the  value  at  the 
wall-layer  edge  is  known.  A plot  of  the  wall-layer  edge  disturbance  is 
useful  therefore  . 

From  (21),  (39),  and  (63), 


1 cH- 

(C/D)(flR)1/2 

r<o  iLl)  , . 

/ (0  + iL)  e e 1 dL 

1 dy 

" 0 2 try  . xll£‘ 

1 e 

J 

“00 

47 


I 


Equation  (10)  then  gives 


U d£ 

r 0y  T)  = 0 


(C/D)(QR)1/2Q 

_2 


2v . RH 
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Y°"  Ye' 

jh:  exp[-  (no  - -hsf)  ] 


Because  the  two  parameters  of  the  right  member  and  YQ  - Y mingle 
frequency  and  position,  they  are  replaced  for  graphic  purposes  with 
parameters  that  do  not.  These  are  Ux^/v,  Uy^/v,  and  u>v/U  (which  is 


equal  to  Q/R).  Then, 


• ■ (?)  ("‘ 


Y0  * Ye  = 


- RTi) 

W)\  v e/ 


Also,  (C/D)(QR)1  ^ is  a function  of  the  R“  (u>v/U2)2  given  in  figure  -4,  and 

2 2 

H and  T]  are  functions  of  R uov/U  given  in  figure  5. 
e e 

In  figure  9,  the  disturbance  source  has  been  located  as  far  upstream 
from  the  leading  edge  as  an  observer  is  downstream  from  the  leading  edge, 
if  his  local  x- Reynolds  number  is  lo"’.  The  location  and  the  frequency 
selected  are  simply  representative.  The  curves  in  figure  9 then  are  those 
that  an  anemometer  should  be  able  to  measure  if  it  were  located  at  the  wall- 
layer  edge  and  translated  along  the  plate.  The  figure  indicates  that,  if  the 
source  is  located  in  line  with  the  plate  (Uy^/v  = 0),  the  amplitude  decays 
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monotonically  downstream.  If  the  source  is  located  below  the  plate,  e.  g.  , 
Uy0/v  = - 1000,  tiie  amplitude  in  the  boundary  layer  above  the  plate  is 
greatly  reduced.  If,  however,  the  source  is  located  on  the  proximal  side  so 
that  its  wake  enters  the  boundary  layer  near  the  stations  observed,  the 
amplitude  peaks  at  a finite  Reynolds  number.  (Generally  two  peaks  should 
be  expected,  corresponding  to  the  two  lobes  of  the  disturbance  in  the  absence 
of  the  plate.  ) The  further  the  plate  is  from  the  source  wake  centerline,  the 
smaller  the  characteristic  disturbance  and  the  greater  the  Reynolds  number 
of  maximum  disturbance. 

In  figure  10,  the  source  has  been  fixed,  and  the  frequency  parameterized. 

4 2 

For  a station  where  Ux/v  s 3 X 10  , the  maximum  disturbance  for  uuv/U 

- 3 2 - 4 _ -x 

10  exceeds  that  for  both  u>v/U  3\10  and  i X 10  . 

In  figure  11,  the  observer  and  the  source  are  both  located  a distance 
10"  v/U  from  the  leading  edge.  The  effect  of  varying  frequency  is  shown, 
lor  Uy q/v  greater  than  1000,  there  is  clearly  a favored  frequency  (if 
perturbation  of  the  boundary  layer  at  that  station  is  desired)  or  a frequency 
range  to  avoid  (if  a quiescent  flow  is  desired).  This  behavior  indicates  the 
possibility  of  "detuning"  a wind  tunnel,  perhaps  by  redesigning  the  turning- 
vane  chord  length  to  avoid  these  particularly  effective  frequencies. 
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Appendix  A:  Limiting  forms  for  (C/D)(fiR) 

The  differential  equation  for  disturbance  stream  function  (20),  can  be 
solved  in  both  the  low-  and  high-frequency  limits.  The  low-frequency 
result,  which  can  be  verified  by  substitution,  is 


[oT|(0)(l  - f. 
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For  evaluation  of  the  ratio  (C/D)  (OR)  , the  behavior  of  o as  'I  — ® is 
needed: 
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where  H (T)  - S)/2.  The  indicated  integration  is  singular  because  G(T]) 
behaves  as 


G(Tl)  ~ (T)  - B)"1  (1  - f_)‘2  4- 


Accordingly,  for  large  T|, 


r ^ 2 

/ fl  - f_(T|  )p  C(T1  ) dT)  ~ log  H + c + 0(11 
Jq  T|  o o o 
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where  c equals  the  finite  part  plus  log  2.  Specifically, 
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- f_(T|  )]2  cen  > dT) 

\\  o o o 


f -pCn  >]2  g(p,  ) - (Ti 
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+ 


Numerically,  c = 2.  125,  Therefore,  at  low  frequencies  near  the  outer  edge 
of  the  boundary  layer, 

[i  + n2  (logH  + c)  + ...1 
2 2 

~ (2Vl )'  1 e”  [1  + cQ2  +...  ] 


The  use  of  (21)  and  (25)  then  gives  the  boundary-layer  amplitude  ratio  as 


(2Yi) (i  - cn2  + 


) 


n « 1.  (A  - 1 ) 


For  high  frequencies,  the  cosh(fi'  ) solution,  (22b),  is  not  uniformly  valid 
because  the  third  term  of  (20)  becomes  important  when  H = (Tj  - 3 )/2  is 
of  order  Q.  Since  the  edge  layer  (Section  3.  3)  adds  viscous  complications 
when  H is  of  order  [in  (TlR)l*  ^ 2,  the  parabolic-cylinder  behavior,  (25),  is 
only  relevant  if  0 < ln(0R).  Equation  (25)  is  essential  for  the  determination 

of  (C/D)(nR!*^,  thus,  this  ratio  will  be  obtained  for  moderately  high 
frequencies  1 « Q2«  fln(w  R>], 
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In  the  notation  of  Abramowitz  and  Stegun  (1966)  when  H is  of  order 
Cl,  o is  a combination  C.yjCH  - 0)  + C Z(T|  - 0)  of  parabolic  cylinder  func- 
tions Vj  and  with  parameter  a = fl"  + 1/2.  The  cosh(flT])  solution  for 
of  order  one  will  match  this  if  is  or(0)exp(03 ) / (2f.)  and  is 
Otj(0  )exp(H0 ) / 2.  It  remains  to  express  yj  and  y^  in  terms  of  standard  solu- 
tions U and  V , so  that  the  asymptotic  behavior  of  y^  and  y 7 can  be  found. 
Since  U « V for  large  arguments,  only  V is  retained. 
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The  use  of  these  and  Sterling's  formula  provide 


1/2  np  ,(1  + 0 )/2  2 2 
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which  with  another  application  of  Sterling's  formula  provides 

C\  /nni1/2  ,rZ/?  x c m + i ,-(^*1/2) 

— ) (CR)  ~ exp  - (U  /2  + wp)  U z 


The  limiting  forms  ( A - 1 ) and  (A-2)  are  plotted  in  figure  4. 


( A-2) 
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Appendix  B:  Edge-layer  solutions 

Four  independent  solutions  were  sought  for  the  edge-layer  equation 


Vhhhh  + ie"h  (°0,  hh  - V = 0 


(B  - 1 ) 


three  of  these,  o*1 1 , 0*2),  and  ©[,3'  were  obtained  by  recognizing 


The  first  inree  or  tncse,  v0  ’ W0  ’ “““  v0 
the  series  expansions  for  large  positive  h. 
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Then,  the  function, 
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is  an  exact  solution,  and 
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Kelvin  functions  in  the  form 
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are  used.  From  these, 
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is  found,  where 
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Similarly,  the  series  representation  of  ©q  is 
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where  Y is  Euler's  constant  and  v(k  + 1)  is  the  psi  (digamma)  function.  Bv 
means  of  a series  manipulation. 
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From  these  equations, 
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is  found,  where 
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Solutions  o{j2)  and  will  be  combined  and  simplified  in  (B  - 36)  and  (B-37). 

(4)3 

An  attempt  to  tind  a fourth  series  Oq  ~ h + • • • was  abandoned  in  favor 
of  asymptotic  representations  of  as  h—  -®.  Equation  (B-l)  should 
possess  two  solutions  that  asymptotically  satisfy  the  "inviscid"  limit,  where 
the  difference  vanishes.  Thus,  solutions  are  sought  of  the  form 
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0^  is  simply  -iOg  , which  precludes  from  being  the  fourth  indepen- 
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dent  solution  desired.  Thus  oQ  is  needed.  It  has  the  asymptotic 
expansion 
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Because  this  series  was  unrecognizable,  it  was  necessary  to  work  with  a 
related  function  <t>  defined  by 


♦ »2<2l»0.ww 


(13-20) 


where  w is  exp  (-h)  here  and  below.  Now  <t>  satisfies  the  second-order 
eq  uation 


w <t>  +■  w<t  +•  iw<t>  2i 

ww  w 


(13-21) 


with  complementary  solution 
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Hq  is  used  rather  than  Yq  because  it  is  directly  expressible  in  Kelvin 
functions : 
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A particular  solution  of  (B-21),  obtained  by  the  vise  of  an  indicial  function  and 
Duhamel  integration,  is 
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By  the  adjustment  of  the  coefficients  A and  B of  the  complementary  func- 
tions each  lower  limit  in  (B-25)  is  set  separately.  Most  convenient  for  * 

y 0 

is  =°.  For  J^,  any  parameter  a is  chosen.  Then,  the  sum  $ + <j>  provides: 
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where  FP  denotes  the  finite  part  of  the  integral.  Two  quadratures  then 
provide  the  fourth  independent  solution: 
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B.  1 Limit  h — °°  (outer  edge) 
(3) 


Asymptotic  forms  for  p^  Pq  *,  and  PQ  are  shown  in  (*3-2),  (B-3), 
and  (B-4).  Use  of  the  expansions 
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results  in  the  following  representation  for 
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whe  re 
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The  two  quadratures  of  <t>/w2  then  result  in  the  asymptotic  form  desired,  i.  e 
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whe  re 
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Although  all  three  e.  can  be  evaluated,  they  have  not  been,  because  they 

are  not  needed  for  zero-order  matching.  The  coefficient  i/3  of  the  cubic 

(2) 

term  in  ( B - 3 2 ) is  needed,  which  justifies  obtaining  on  and  its  limiting 


form  shown. 


B.  2 Limit  h — -<*>  (inner  edge) 


F rom  (B-5)  it  is  clear  that  as  h — -<*>,  Oq  ~ i exp(-h),  and,  from  (B-18) 
that  U)oQ  ~ exp(h).  Examination  of  (13-11)  indicates  that  Oq2'  grows  as 
h -»  -cd  much  more  rapidly  than  exp(-h). 

The  same  is  true  for  There  is,  however,  a linear  combination  of 

Oq*- ' and  o<3)  for  which  this  rapid  variation  cancels.  That  combination  is 
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If  x Ze  f can  be  written  as 


whe  re 
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The  integral  in  (B-37)  vanishes  as  x — 0,  since  kr(x)  - iki(x)  is  of  order 
x In  x.  Equation  (B-37)  does  not  contain  any  of  the  (xe  ) terms  that 

caused  and  to  become  unbounded  for  large  x.  Since  ' (xe^^) 

vanishes  for  large  x (the  inner  edge),  an  asymptote  for  f that  is  merely 
exponental  in  h is  left. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer. reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  n polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  • emistry.  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materia's,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics,  quantum  electronics, 
lasers,  and  electro-optics,  communication  sciences,  applied  electronics,  semi- 
conducting. superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment,  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetosphe nc  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere ; solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy,  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere,  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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